Many-body dynamics of Rydberg excitation using the i7-expansion 
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We investigate the excitation dynamics of Rydberg atoms in ultracold atomic samples by ex- 
panding the excitation probability and the correlation function between excited atoms in powers 
of the isolated atom Rabi frequency fi. In the Heisenberg picture, we give recurrence relations to 
calculate any order of the expansions, which ere expected to be well-behaved for arbitrarily strong 
interactions. For homogeneous large samples, we give the explicit form of the expansions, up to Q!^, 
averaged over all possible random spatial distributions of atoms, for the most important cases of 
excitation pulses and interactions. 

PACS numbers: 32.80.Rm, 03.67.Lx, 34.20.Cf 



I. INTRODUCTION 

In recent years, the interactions between ultracold Ry- 
dberg atoms have been extensively studied for their possi- 
ble significance in molecular and ultracold plasma physic, 
as well as quantum computing. It has been proposed to 
use strong interactions between Rydberg atoms to entan- 
gle neutral atoms and achieve fast quantum gates [l|, Q ■ 
One effect of these interactions is the excitation block- 
ade, where one excited atom prevents the excitation of 
nearby atoms. This effect could be utilized to realize 
scalable quantum gates In macroscopic samples, the 
blockade mechanism manifests itself as a suppression of 
Rydberg excitation. A local blockade of Rydberg state 
excitation in a mesoscopic sample due to strong van der 
Waals (vdW) interactions has been observed [4j] using a 
pulse-amplified single-mode laser. We have proposed a 
mean-field type model to explain these experimental re- 
sults [4] . In the model, a distribution of mean-field shifts 
was constructed for which a distribution of excitation 
probabilities was calculated. The agreement between the 
theoretical model and experimental measurements was 
good. Suppression of Rydberg excitation has also been 
measured in two-step excitation processes using cw ex- 
citations 0, B 0] ■ The atom counting statistics of Ry- 
dberg excitation is significantly modified by the interac- 
tions and its sub-Poissonian character has been observed 

Strong Rydberg-Rydberg interactions can impose cor- 
relation between atoms within a range of a few /im so that 
many-body treatments are, in principle, necessary. Vari- 
ous treatments [1, [§, [13] have been proposed to describe 
these systems. For example, rate equations were used for 
certain two-step excitation schemes in [1, [ll|. However, 
this approach is not applicable for the systems we pri- 
marily consider because spontaneous decay and related 
decoherence effects are negligible in our case. A differ- 
ent approach is to numerically calculate the many-body 
wave function [l^, [13] ■ For strong correlations between 
nearby atoms, the concept of pseudo(super)atoms can 
be introduced (lo| to reduce the number of many-body 
states needed for the numerical simulations. We can nu- 
merically evaluate the many-body wave function of meso- 



scopic systems of ~10 /im diameter [l2|. Large systems 
are difficult to describe in this method, so that the in- 
teractions with faraway atoms is modeled by mean-field 
shifts while the correlations between nearby superatoms 
are fully accounted for [l^]. 

In this paper, we show how to calculate the f2- 
expansions of the excitation probability and correlation 
function in the Heisenberg picture. The explicit forms of 
the lowest orders of these expansions are derived and for 
large systems their ensemble averaged counterparts are 
given. 



II. n-EXPANSION IN THE HEISENBERG 
PICTURE 

We consider a system of N ultracold two-level atoms 
for which the upper level is a Rydberg state. We as- 
sume that the thermal motion of the atoms is greatly 
reduced and can be completely ignored, leading to the 
so called "frozen" gas approximation. The interactions 
between Rydberg atoms are presumably strong so that 
many-body effects may occur. 

We start with the many-body Hamiltonian of interact- 
ing two-level atoms (with h = \) 



N 



N 



H = A^aL + -^HOa:,+u;*(i)<Tg 



N 



i—l,j>i 



(1) 



where A is the frequency detuning from resonance and 
Kij are the interaction strengths between Rydberg atoms. 
The second term in the Hamiltonian is the dipole oper- 
ator representing the interaction with the optical field. 
The function w{t) is the time evolution (envelope) of 
the laser pulse. The cr-operators are defined as ct^^ = 

\o!i) {Pil (Xi P^^~^\ where a, (3 refer either to the ground 
state g or the excited state e and is the identity 

operator in the subspace which is an orthogonal comple- 
ment to \gi) {g^ \ ® \ei) {ei\. 
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Solving this Hamiltonian is trivial in two limits, when 
all Kij —^ oo and when all —^ 0. For Kij —f and 
real w{t), we get the isolated atom excitation probabil- 
ity given by the Rabi formula Poxc(i) = sm'^{il.W{t)/2), 
where ilW{t) — J^^ w{t')dt' is the pulse area and to 



is the initial time of laser excitation. For k.; 



oo, 



an arbitrary big ensemble of atoms becomes fully block- 
aded. This means that there cannot be more than one 
excited atom in the sample. Such systems are effec- 
tively two-level systems, with the collective ground state 
\G) — 13152 • ■ -Qn) and the collective excited state \E) — 
1/Vn \gi . . . Ci . . . (?Ar). Consequently, the equivalent 
Hamiltonian H' in the limit kh — j-cx) is 



Act 



Nil 



EE 



{w{t)aEG+w*{t)aGE)- (2) 



This is formally the Hamiltonian for isolated atoms, only 
the Rabi frequency scales as v^- If the detuning A is 
zero and w{t) real, the solution is 



1 



Nnw{t)/2 



(3) 



We can solve the equations of motion by expanding the 
CT-operators in powers of fi. It follows from Eq. ([3|) that 
the expansion is well defined and convergent in the limit 
Clearly, since it is also fine for isolated atoms 
(the limit Kij—fO), we expect that it is well defined and 
convergent for arbitrarily strong interactions. 

It is remarkable that, in principle, all terms in the 
expansion can be calculated exactly. In this sense, the 
Hamiltonian ([1]) is exactly solvable. The usual problem 
with a many-body Hamiltonian is that the nonlinear part 
containing interactions is very difficult to solve. Here 
the nonfinear part Hq = AE^^ct^^ + ct^^ct^^ 
is exactly solvable. This Hamiltonian cannot change the 
number of excited atoms and thus any pure collective 
state with a fixed number of excited atoms is an eigen- 
state of Hq. This is essentially the reason allowing the 
calculation of any term in the fJ-expansion. 



A. fi-Expansion of the excitation probability 

The evolution of the CT-operators in the Heisenberg pic- 
ture, governed by the Hamiltonian ([T]), is given by the 
following equations 



[9it)Kg - 9* {t)^le] ^ 



dr ^ 2 
dal 



(4) 



eg 



dr 



Aalg + z-g*it) [2CTL-l]-fz5]^.,CT:,CTL. (5) 



These equations can be simplified by removing the A 
term using new scaled (dimensionless) variables 



T = t/T, LO = QT, ai^e 



3 p-i^t 



6 = AT, fcy = KijT, /(r) = g{TT)e 



eg' 

iSr 



where T is the puls e dura tion. For Gaussian pulses, we 
use T = Tpyy HM /V2\ii 2, where TpwHM is the interval 
between the instants for which the laser power is a half 
of the maximum value. The equations for the new CT- 
operators can then be rewritten as 



dal 



dr 



= ^[/(rK,-r(r)CTy , 



(6) 



da, 



eg _ 
~"2 



.rw[2CTL-i]+z^%CT:^CTL. (7) 



We solve the equations of motion by expanding the 
CT-operators in power of uj. 



^ ee ee ' ee ' ^ ^ ee ' 

(7 — fJ + UJ(7 -\- to (7 + 



(8) 
(9) 



The w-expansion of {al^ir)) ((CTeg(r))) contains only 
even (odd) power of w if all the atoms are initially in their 
ground state. This can be shown as follows. Instead of 
\g) and |e) we can use \g') and |e') defined as 



\9) 



(10) 



where the constant phase factors satisfy fe — ^g — i"". 
We can then define new operators CTe'e' and a^'g'- It 
follows from the definitions of these operators that for 
any r and uj the following relations have to be satisfied 

ae'e'{T,Uj) = aee{T,U!), CTg'g' (t, w) = -CTeg (t, tj) , (11) 

which gives 



{ae'e'{T,Uj)) = (CTee(T,w)), 
(ae'g'{T,Uj)) = -(CTeg(r,w)). 



(12) 



Similarly to Eqs. ([6]) and ([7]), we can write evolution 
equations for the new operators CTe'e' and ag'g'- We can 
also write Eqs. (O and ([7]) for the opposite sign of omega 
Lu — !■ ~uj. However, after the substitution lu' ~ —uj, the 
time evolution of the expectation values {aeeit,uj)) and 
(CTeg(t, w)), and the time evolution of {ae'e' it,uj')) and 
{(Je'g' {t, uj')) are given by the same differential equations. 
Since the initial conditions are the same, the solutions 
have to be the same as well 



{aee{T,Uj)) = {ae'e'{T,Uj')) = {ae'e'{T,-Uj)), 
{aegiT,Uj)) = {ae'g'iT,Uj'))^{ae'g'{T,-Uj)). 

Combining Eq. p^ and Eq. p^ we conclude that 

(CTee(r,w)) = (CTee(T, -w)), 
(CTeg(T,U;)) = -(CTeg(r,-w)). 



(13) 



(14) 



Therefore, if all the atoms are initially in the ground 
state, (CTee(i, w)) is an even function of uj and {aeg(t,uj)) 
is odd. All these symmetry properties can be checked 
against the results below. We note that this statement 
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is not true for arbitrary initial conditions or for the a- 
operators themselves. 

According to Eqs. ^ and ([7]), the differential equa- 
tions for (Tee"'' and CTeg"^ are 



da, 



dT 



(15) 



da, 



i{n) 
eg 



.nt) 



n-l 



dT 



j^i p=0 



ee 7 



(16) 



where (5„i is the Kronecker delta function. To start the 
recurrence procedure, we need the initial operators al^e^ 
and al^gK They are derived from Eqs. ([6]) and ^ for 



<^;r(T)=<^r(ro), 
ai^^\r)^alf\To)Q\T~To), 

where tq is the initial time and is defined as 



(17) 



Q\t) = exp i r'^ki 



(18) 



The equivalent integral form of Eq. (I15p is 

/(n)4"-^nn) 

-r(ri)<7^i"-i)(ri)". 



Ki^'Hr)^- I dT, 



(19) 



We get an equivalent integral form of Eq. pB]) in two 
steps. After multiplying this equation by (5*(— t), group- 
ing terms with iJeg"' and utilizing the time independence 
of a ee*^^ , another differential equation is obtained 



«/*(n) 



(2a:("-i) - 5„i) 



Q\t). (20) 



The corresponding integral equation is 



4"Ht) = / dn 



To 



^— (2CTe^e '[Tl)~6nl) 



n-l 



\t-t,). (21) 



+ *E%E4'H^l)'^ei"-^^(^l) 



In our approach the recurrence relations (|19p and ([2T|) 
provide the way of calculating any order of the expan- 
sions. To calculate (Teg"\ one needs to calculate CTei"^ 
first. What we really want is the expectation values of 
the CT-operators, especially {a\^). With the help of the 

recurrence relations, we can express all a^^e^ and iJeg"^ 
in terms of tr^i"^ and at^g^ . Since we can easily find the 



expectation values of any product of (T^e°^ and aig'' , in 
principle, we can find exactly any term in the expansions 
([8|) and ((9]). In our calculation, we assume that all atoms 
are initially in the ground state. Consequently, the fol- 
lowing expectation values are necessarily equal to zero 
for any atom i 

{alf. ..) = (... ap ) = (a]^\ ..) = (... a^l") ) = 0. (22) 



=.i(o) 



The operator Q only produces phase factors. For ex- 
ample, for any a, the following relation are satisfied 



dQ\T)) = (a) , 

aQiT)a^e^"^) =exp(zT%)(aa^f 



(23) 



They follow from the definitions of Q and the a- 
operators. Will use this property of Q to have all kij in 
our formula exclusively in the phase factors exp(irA;ij). 

Since in the Heisenberg picture the wave function is 
time- independent, the average ( ) can go through inte- 
grals contained in the recurrence relations. We find a 
new recurrence relation for ^(Tel"'^ after substituting the 

~i{n-l) 
Teg 



expression ([2T|) for a. 



into Eq. H 



Re 



JdT.^^^^l^^ (r(n)(5„.-2(a:("-)(n))) -Y.k,j:{a%^\n)ai^-'Hri 



(24) 



j^i p=0 



where F{t) = JJ^ dT' /{t'). For n = 2, this equation, together with expressions (|22p . leads to 

r r 

[^^eHr))^ J dri J dr2<i>(ri,T2), 



(25) 
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where 



For any symmetric fmiction ^'(Ti,r2) — ^'(T2,ri), the 
foUowing relation is true (if the integrals exist) 

T ri T T 

J dTl j dT2-^{Tl,T2) ^ j dTl J dT2*(Tl,T2) 

To To To Ti 

T r 

= ^ y dri y dT2^'(ri,T2) . 



Using this property, we finally get 



(27) 

The last relation shows that the effects of interactions 
comes through higher orders than n — 2. This, of course, 
would not be true if the initial conditions were different; 
in our case, since there are no excited atoms initially, 
there are no effects due to Rydberg-Rydberg interactions. 
For these initial conditions, the excitation always starts 
as isolated atom excitation. 

Besides Eq. (f24| . we can derive more auxiliary rela- 
tions for expectation values. They are useful if the same 
type of recurrence relations is repeatedly used. In our 

case, the expectation value (^S'(r')o-ei^''(r)^, for an arbi- 
trary operator S{t'), will be particularly helpful. After 
substituting expression for n ~ 2, and then using 
Eq. pT|) for n = 1, we derive 



\F{r) 



S{r') 



^!^(5(r')a:W(ro)a^f (to) 



X / dTie'(^i-^«)'=-F(Ti)(i^(T)-F(Ti)). (28) 



The complex conjugate of this formula leads to the ex- 
pectation of ((Tee^^(T)S'(r') ) for an arbitrary S{t') 



X / dTie-'(^i~^°)'''^F*(Ti)(i^*(T) - F*(ti)) . (29) 



Calculating ( al^e^ (r) ) takes more effort, so we split its 



contributions, given by Eq. (j24p . into three parts 

[Ki'\r))^^h,-h2^h,. (30) 



The integrals are defined as follows 



/41 = / dri^- — ^-^-Lrc 



'42 



nn){FiT)^Fin)) 
F{t)-F{t,) 



(31) 



lis. = ^ fey Re 



To 



dTl 

2 

F(r)-f(ri) 
dn ^ 



(32) 



(33) 



The integral /41 can be calculated easily since ( al^g'^ (r) 



is given by Eq. ((27)) . After integration over T2, we obtain 



/41 = 



\nr)\' 

16 



- Rc 



Fir) 



T 

JdTig*{Tl)F\Tl) 



(34) 



Only /42 and include interactions. The calculation 
procedure is to apply Eqs. p9)l and pT|) repeatedly, until 



the only operators left are al^g"^ and , whose expec- 
tation values are trivial to calculate. Using the auxiliary 
relations ([24|) - (|28p often simplifies the derivation of ex- 
pectation values. To evaluate \(Tl^g\Ti)aii^\Ti) \ in 142, 



(35) 



we can use Eq. ([28)) . The result is 



J42 =^fcylm 



F{T)Jdn ^ F (ri 



This integral is canceled out by one of the I43 terms, so 
there is no need to consider it in detail. 

It is convenient to have kij in the exponential (phase) 
factors only. This can be done for the sum X4 = 142 -I- 1^3 . 
After several partial integrations, the simplified form of 
T4 is 

- r 

^4 =iEi^<^ / ^^^1/(^1) inr)-mri)) 



dr. 



2g*(T2)F*(T2)(e^(^i-^^)'=--l) 



(36) 



For a system of N atoms, one just needs to calculate 
these integrals for given kij . 

The derived expressions for /41 and Z4 are formally 
sufficient but they might not be very convenient in the 
limit of a large number of atoms N , where we replace 
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the sum ^j^i by the integral J pd^R, where p is the 
atom density. For the terms we exphcitely consider in our 
expansion, it can be shown by direct calculation that this 
replacement is equivalent to averaging over all possible 
spatial distribution of atoms. For a large sample and 
an arbitrary pulse shape, an easier way is to first find 
the sum in Eq. ([55)1 (i.e. the integral J pd^R). The 
last form of 1^2 + ^43 is very convenient to account for 
any angular dependence of /cy . As an example, we give 
the result for the excitation probabilities in large systems 
(surface effects are ignored) with uniform densities. For 
large homogeneous systems and kij = 2-kCsT / R'^ , where 
Cs may be angular-dependent, we obtain 



1 =A(|C.|T(ri-T2)) 



3/s 



In this expression, A is a parameter which depends on the 
type of interactions. Note that n — T2 is always positive 
in X4 . 

The ensemble averaged X4 for large homogeneous sys- 
tems is 



J4 =A p {CsTf/'^c 



Jdnf{n){F{T)-2F{n)) 

To 

jdT2g*{T2)F*{T2){Ti-T2f/' 



(37) 



For resonant excitation and real /(<), the ensemble aver- 
aged expansion of excitation probabilities is 



(38) 



where /gat is the saturation laser intensity (for iso- 
lated atoms) and T is the pulse duration. The intro- 
ductio n of / gat allows using a single formula Pexc = 
s\T? {^J I / IsatT^ / 2) for rcsouant excitation of isolated 
atoms regardless of the type of excitation pulses. For 
s = 6 we have the van der Waals interactions and for 
s = 3 the dipole-dipole interactions. The values of the 
parameter 7 in various cases of laser pulses and interac- 
tions are presented in Table |T1 We note that the parame- 
ter 7 for the angular dependent dipole-dipole interactions 
(with aligned dipole moments) V{R) — — cos^ ff) is 

7 — 47'/3-\/3, where 7' corresponds to the isotropic in- 
teraction Uz/R^- 

These formulae have the first contributions of inter- 
actions to excitation probabilities. Note that we never 
assumed that the interactions were weak. Actually, Eq. 
P6p is consistent with the limit kij — > c» for any i, j. In 
this limit, there is no contribution from the exponential 
terms since they oscillate infinitely fast. The remaining 
part reproduces exactly what one gets when the exact 
solution ^ in the limit kij — s- 00 is expanded in powers 

of Ld. 



TABLE L Tlie parameter 7 in tlie expansion (I38|l of excitation 
probabilities for various interaction potentials and excitation 
pulses. The pulse envelope for a Gaussian pulse is gij) = e~'^ 
and for a square pulse is gij) — 0(1 — r). 



pulse type 


7 




^(l-3cos2 6) 




Gauss, pulse 


32.1138 


24.7212 


10.8627 


Square pulse 


5 


15%/3 


1287r^ 
189 



B. f2-Expansion of the correlation function 

We also calculate the spatial correlation function 
P{i,j) between different atoms i and j 



{Ke{r)aie{r)) 



(39) 



We first expand c(i,j) = {(^\ei'^)^ie{'^)) powers of 
u). Since {<j],e{'^)'^ee{''')) is essentially a probability, we 
expect only even terms in the expansion 

{Ke{r)aL{r)) =u:\^^\^,J)+u:^c^^\^,3) + ... , (40) 

where c'^'^''(«, j) = due to the initial conditions ^F^. We 
can directly verify that the lowest odd c'^^\i,j) terms 
vanish. According to Eq. ([^ . c^^^(«, j) is equal to zero 



c«(z,j) = (a:Wa^(°)) + (a:(°)a^W)=0. (41) 



Using Eq. (I19|) . we also find 



\F{t 



— "i,— 0, (42) 



since i ^ j- Relation ([27ll also follows from the last for- 
mula because = ^leO'le- We can exphcitly show that 
c(3)(i,j) = (^ali'\T)aii'\T)) + {ali'\T)aii'\T)) ^ 

using Eqs. 1^ and ^ with S = alf^'' and (^alf^^^ = 

0. Therefore, c'^^^(i,j) is the first nontrivial term in the 
expansion of {al^iT)ai^{T)) 



{^,J)= (a^W {r)al^ (r)) + (af ) {r)al^ (r) ) (43) 



This is expected since the initial excitation probability of 
any atom is ~ w^. The first and last term in the right- 
hand side of the last equation are very similar. If we 
know one of them, it is easy to find the other one. To find 
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alf\T)aif\T)Y we apply Eq. ^ for S{t) = a2\T) 



16 



X / rfrie*(^i^^'')'=^=i^(Ti)(F(r) - F(ti)) . (44) 



The expectation value in the last formula is easily evalu- 
ated using Eq. together with the following identity 



■ ge yv/"ge V'uy"eg V'u;"eg V'uy/ "is-'jq, 

where we assume i ^ j. The result is therefore 



16 



16 



T 



(45) 



A more convenient form obtained after partial integration 



|F(7 



16 



1 

+ 16 



T 

j dne'^'''^^F{n){F{T)- 2F(ti)) 



(46) 



Instead of using Eqs. (fTO|) - (PT|l . another approach to 
find \d-e^e\T)dii'^\T)\ is to utilize the identity a^^ = 



^L^ie^ which gives aif'' = alf^aif^ + a'ee'' cr'Jk^ + 
di^^ ai^e'' + ^ee'' (^ee'' ■ The expectation values of the first 
and the last term in this expansion vanish so that 

+ {K^e\r)o'}eHr)<^i'^e\r)). (47) 

To evaluate (a^ee^ {T)(Ti^e^ {T)(T'ee^ (t)'^: wc first apply the 
recurrence relation (fT51) on o'ee^'' (t) anddS"(T) 

r 

X — J dTie-'(^i-^°)'='V*(Ti). (48) 

To 

To get the expectation value in the right-hand side of the 
last equation, we use Eq. ([^5]) together with the following 
identity 

(4°^(ro)a^i°)(ro)a^(°Hro)a,f )(ro)) = 6... 



The result is 



To 

dT2 e*("^-^i''=-F(T2)(^^(T) - F(t2)). (49) 



This can be additionally transformed using partial inte- 
gration to get 



-F*{t) 
16 



dTl/*(Tl) 



dT-2e^(-^— i)'=-F(r2)(F(T) - 2F{t2)). (50) 



The only remaining term to calculate is the last term in 
Eq. (|47| . From Eq. (fT9)l . we immediately find 



\nr) 



^i^T(jo)a^{r)a^(r^. (51) 



We proceed using Eq. to express (7^i^''(T) in terms 
of o'eg^'' and Ug'^t^ , and then apply Eq. pT|) on (Teg^'' and 
(T^i^^. The simplified form after these substitutions is 



-^\F{r) 
16 



X % dT,r{T,) / dT2e''(^^-^i)'=-F(r2). (52) 



Partial integration will leave fcy only in the phase factor 



\F{r)\- 



\Fir)\^ 
16 



r 

Jdnr{n)Fin) 



+ -^^ d^in^i) /dT2e'(— ^)'=-i^(r2). (53) 



16 



Combining Eqs. P5)l . ([50)1 and we get our final 

formula 



Jdne'^'''^^f{Ti)F{n) 



(54) 



The first term in the expansion of the product 
{^lei^}) i^eei''')) also proportioual to uj^. According 
to Eq. (gll), it is 



16 



(55) 
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From Eqs. ((39|) . (|54p and ((55)) . we find the correlation 
function for low laser power to be 



\F{r) 



(56) 



This simple formula includes the effects of interactions, 
frequency detuning, and possible frequency chirp. As 
pointed out in [l^, [I4I , the numerical calculation of the 
correlation function is much more demanding than the 
calculation of the excitation probability in the super- 
atom approach. Both calculations include running sim- 
ulations for many random spatial distributions of atoms 
to find ensemble averaged probabilities and correlation 
functions. In the superatom approach, the conceptual 
issue in the calculation of the correlation function is re- 
lated to the fact that superatoms are extended objects. 
The question is how much the numerical correlation func- 
tion averaged over many random spatial distributions of 
atoms is accurate over distances less that the average size 
of a superatom. 



III. RESULTS AND DISCUSSION 

We know from theoretical 0, and experimental 
studies [1, that excitation fractions are easily sat- 
urated in large systems of strongly interacting atoms. 
This, of course, cannot be concluded from the first two 
terms in Eq. ((38)) . However, assuming that the prob- 
abilities are saturated, we can use the maximum exci- 
tation probability, determined from these two terms to 



estimate the saturated excitation fraction. From the cor- 
responding laser irradiance, we can get an estimate of the 
saturation intensity /q for which the excitation fraction 
becomes saturated due to interactions. 

We introduce a suprcssion factor Nd defined as the ra- 
tio of the maximum excitation probabilities for isolated 
atoms and for interacting atoms. This Nd is basically 
the number of atoms in the largest possible region con- 
taining no more than one excited atom. The maximum 
excitation probability Pq obtained directly from the first 
two terms in Eq. p8|) is 



3/4 



1 + 7P(|C,|T)3A 



(57) 



Obviously, this formula underestimates the excitation 
probability in the absence of interactions Cs — > since 
it gives Po — 3/4 instead of Pq = 1. This is due to the 
truncation of the expansion. Also, in the limit Cg — > 00, 
the first two terms in the expansion of Pexc, according 
to Eq. (121), lead to the estimate Pq = 3/4 A'^ instead of 
Pq = 1/A^. It seems that, in the first approximation, the 
limitations due to using truncated expansions could be 
overcome by replacing the factor 3/4 by 1. This replace- 
ment leaves the estimated excitation fraction, as a func- 
tion of the interaction parameters, basically the same. 
Another procedure, based on Eq. ([3]), leads to the same 
correction. We can use a modification of Eq. ([3]) to get a 
model dependence which is correct in both limits of Cs ■ 
The simplest modification is 



Po. 



sm 



./NdI/IsatTT/2) /Nd 



(58) 



This Nd is equal to one for isolated atoms and Nd = N 
for a fully excitation blockade. For arbitrary interactions. 




FIG. 1: Rydberg excitation fraction as a function of laser 
intensity I. The solid line is obtained using Eq. (|59p for 
I > lo and tiie model dependence (|58p for I < Iq. The result 
is comparable with the numerical solution (dashed line) of 
the many-body wave function [l^. The dotted line is the 
noninteracting limit. The parameters used for this simulation 
are Ce = 2.64 x 10^^ x 7/60, p = 6.5 10^" cm~^ and a laser 
bandwidth T = 120 MHz. 
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FIG. 2: Density dependence of the Rydberg excitation frac- 
tion. The solid line is this dependence obtained using Eq. 
(|59p and the dashed line is obtained using a mean field model 
described in 0] for ///sat ~ 0.2. The values of Ce and F in 
this figure and Fig. [1] are the same (the pulse duration is 
adjusted to get this bandwidth). The highest atom density 
Pmax in this figure corresponds to the atom density in Fig. m 
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FIG. 3: Correlation function between excited atoms for var- 
ious detunings and chirp parameters. In (a), the laser band- 
width r was varied as indicated, assuming that the excess 
bandwidth is caused by linear chirp. The solid line (F = 60 
MHz) represents the dependence with no chirp. The dashed 
line crossing the solid one is the only dependence with a neg- 
ative chirp. Characteristic particle correlations exist only for 
a positive chirp. In (b), the frequency detuning was varied 
as indicated, assuming no chirp (F = 60 MHz). Since the 
interactions are attractive, characteristic particle correlations 
exist only for positive detunings. This figure demonstates the 
strong dependence of the correlation function on detuning. 



is obtained requiring that the first two terms in the 
expansion of Pcxc given by Eq. (I58p match the evalu- 
ated terms in Eq. ((38|) . Strictly speaking, this model 
dependence can be only used for low laser power, or 
more precisely, up to the first maximum of the excita- 
tion probability. Due to the saturation nature of the ex- 
citation probability, we assume that this maximum is a 
reasonable estimate of Pq ■, even for somewhat higher laser 
power. The excitation fraction obtained using these as- 
sumptions agrees well with the numerical solution [lol | 
of the many-body wave function, as shown in Fig. [T] 
Therefore, the estimates of the saturated fraction Pq and 
the saturation intensity /q for interacting atoms are 



Pn = 



^sat 



1 



1+7P(|C.||T)3A 



(59) 



This formula can be also viewed as a density dependence 
of the saturated excitation fraction. This dependence 
is shown in Fig. [2] for Ce = 2.64 x x 7/60 and a 
laser bandwidth F = 120 MHz. The largest atom density 



Pmax in Fig. [2]is 6.5 x lO"'^" cm~^. In this figure, we also 
show the result from a mean-field model described in [3] ■ 
The laser power used in the mean-field calculation was 
I /Isa,t = 0.2. The agreement is fairly good. For the the 
above values of Cg and F, and p — Pmax, we also have 
the result on saturated excitation fraction from many- 
body simulations fio'l . The saturated excitation fraction 
from the full numerical calculations slowly varies from 3.7 
to 4.0 %, while we get Pq = 3.7 % using Eq. The 
same parameters were used to get the Rydberg excitation 
fraction in Fig. [TJ The only difference is that we adjust 
the pulse width to get the experimental F — 120 MHz, 
while in [j] and [10| this bandwidth was caused by a linear 
chirp. We can compare the prediction of our formula ([55]) 
with the results from the numerical treatment for the 
experimental parameters in [5]. The simulation Jjy was 
done for Cq = 4.97 x 10^^, T = 37.5 ns and a sample 
density of 2 x 10® cm^'^. The resonant excitation fraction 
was about 7.5 %. For the same parameters, we get 8.2 
% using Eq. ([5^ . It is interesting that Eq. ([5^ gives 
results very similar to the results of the full numerical 
treatment in both simulations, even though the densities 
and pulse durations are very different. 

In Fig. Ela) we show the correlation function obtained 
from Eq. ([55]) for various laser bandwidths F caused by 
a linear frequency chirp. The pulse width of a Gaussian 
pulse is 7.3 ns and the interaction is characterized by 
Ce = 2.64 X 10^^ x 7/60. the same parameters used in 
Q and the simulation [l^. The dependence for F = 120 
MHz can be compared with the numerically obtained de- 
pendence in [l3] for the lowest optical field amplitude. 
These two dependences are remarkably similar. This fig- 
ure confirms that there is a region where the correlation 
function is greater than one only for a positive chirp (the 
laser frequency linearly increases during excitation). In 
Fig. Eljb) we show how the correlation function is affected 
by frequency detuning. This figure demonstrates that the 
correlation function is very sensitive to frequency detun- 
ing. Although the detunings used in this figure are just a 
fraction of the laser bandwidth, the correlation function 
dependence varies significantly. Positive detunings par- 
tially compensate the effect of attractive interactions so 
that certain separations between excited pairs are prefer- 
able for a given positive detuning. Therefore, the corre- 
lation function is greater than one at such separations. 
This is not the case for negative detunings. These prop- 
erties are also shown in numerical simulations [To| . From 
these figures one may expect that relative small negative 
detunings can easily cancel out the effect of a possible 
positive chirp, which is easy to show using Eq. (j56p . It 
is worth noting that the relatively simple formula (j56p 
reproduces well the main results of demanding numerical 
calculations. 

Even though our formula for excitation probabilities 
([59)1 and the correlation function ((56)) are primarily de- 
rived for lower laser power, they describe related physical 
phenomena reasonably well despite their simplicity. Cou- 
pled with the fact that numerical many-body calculations 
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are very demanding, these relatively simple formulae can 
be useful for actual experiments. 

IV. CONCLUSION 

We have investigated the many-body excitation dy- 
namics in ultracold Rydberg systems using the Q- 
expansion. We have shown that the equations of motion 
can be solved by expanding the (T-operators in powers 
of r2. Different terms in the power expansions of the 
cr-operators can be calculated using the recurrence rela- 
tions. These recurrence relations are used to evaluate the 
expansion of the excitation probability up to the term. 
For homogeneous large samples, the expansions obtained 



are additionally ensemble averaged in various cases of ex- 
citation pulses and interactions. We have also derived an 
explicit form of the correlation function between excited 
atoms for laser power and compared it with the recent 
numerical calculations. Our derived formulae agree well 
with recent numerical many-body simulation [lO| . Their 
simplicity makes them easily applicable in a range of ex- 
perimental conditions 
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